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1 

Method and apparatus for backscatter absorption spectroscopy 

The present invention relates to a device for the determination of a physical 
S feature of a medium, comprising: 

- a light source for producing radiation; 

- a probe with at least a first and a second optical fiber, the first optical fiber having a 
first diameter and being arranged to deliver the radiation on a sample of the medium 
and to collect first backscatteied radiation firom the sample, the second optical fiber 

10 having a second diameter and being arranged to collect second backscattered radiation, 
wherein the second optical fiber is positioned alongside Ihe first optical fibei^ 

- a spectrometer for producing a first signal based on the first backscattered radiation, 
and for producing a second signal based on tiie second backscattered radiation; 

- a processor arranged to determine a differential backscatter signal usmg the first and 
IS second signals. 

Such a device is known from Amelink et al [1]. There, the device is used to 
determine particle sizes in superficial layers. The device is suitable for measuring 
particle sizes in for example an aqueous suspension with polystyrene spheres, but is not 
fitted to accurately measure particle sizes in Uving tissue. So, determining whether 
20 living tissue is normal or precancerous, by way of measuring particle sizes in living 
tissue is not veiy promising. 

« 

In Doombos et al [2] the optical properties of human tissue are determined in 
vivo using a spectroscopic arrangement with ten optical fibers. One of the fibers is used 
to irradiate a sample, and nine other fibers collect the reflected light. By itsing a 
25 multitude of fibers to collect the reflected light, it is possible to calculate scattering and 
absorption coefficients of the sample. However, the.method is not suitable for locally 
measuring the optical properties of tiie tissue. In particular, only mean values of the 
absorption coefficient of a relatively large part of the sample can be .detenmned. 

It is an object of the present invention to locally measure a physical feature, such 
30 as a concentration, of a substance in a medium. 

The object is achieved by a device as described above, characterized in that the 
processor is arranged to calculate the physical feature using a predefined mathematical 
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model, the differential backscatter signal and a curve fitting mechanism. Witti the 
invention very local measurements of for example concentrations are possible. 

In an embodiment^ the first and second diameter are selected depending on a 
mean free path of photoiis sent into the sample. Very good agreement between 
5 measurements and two mathematical models are obtained when relating the fiber 
diameters to the mean j&ee path of photons sent into the sample. 

The fiber diameters may differ, but preferable they are equal and referred to as 

selected fiber diameter. 

According to an embodiment of the invention, in the predefined model mentioned 
10 above, the differential backscatter signal is an exponential function of two times the 
mean free path, and wherein the selected diameter is selected so that: 

where djiber is the selected diameter and mfp is the mean free path. 

In another embodiment, in the predefined model mentioned above, the 
1 5 differential backscatter signal is a ftinction of the selected fiber diameter, and wherein 
the selected diameter is selected so that: 

where dftb^ is the selected diamet^ and mfp is the mean firee patii. 

In a special embodiment, the medium as mentioned above is living tissue. Since 
20 the device can be used to very locally measure concentmtions, it is very suitable to 
detect (pre-)cancerous tissue by measuring excessive oxygen consumption 
accompanied by a low capiUary oxygen saturation. 

The invention also relates to a method for the determination of a physical feature 

of a medium, comprising: 
25 - producing radiation with a light source; 

- placing a probe with at least a first and a second optical fiber, the first optical fiber 
having a first diameter, on a sample of the medium; 

- placing a second optical fiber having a second diameter, alongside the first optical 
fiber; 

30 - sending light coming from the light source, through the first optical fiber; 

- collecting first backscattered radiation through the first optical fiber and second 
backscattered radiation through the second optical fiber. 
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" producing a first signal based on the first backscattered radiation^ and a second signal 
based on the second backscattered radiation; 

- determine a differential backscatter signal using the first and second signals, 
characterized by calculating the physical feature using a predefined model, fbo 

5 differential backscatter signal and a curve fitting mechanism. 

Preferably the first and second diameter are selected depending on a mean firee 
path of photons sent into the sample. 
The method may also comprise: 

- simultaneously measuring backscatter radiation on different locations of the sample ; 
10 - determining a physical feature for the different locations; 

- calculating a standard deviation of the physical feature. 

By doing several local measuremrats» heterogeneities oj^ for example, 
concentrations of substances in the sample can be discovered: 

Furthermore, the invention relates to a computer program according to claim 16. 
15 Finally, the invention relates to a data carrier according to claim 17. 

The present invention will be described below with reference to exemplary 
embodiments and the accompanying schematic drawings, in which: 

20 Fig. 1 is a schematic diagram of a measuring device according to a preferred 
embodiment; 

Fig. 2a and 2b show cross sections of a sample and two fiber tips in the situation 
wherein the mean firee path of the photons is much larger than the diameter of the 
fibers; 

25 Fig. 3 shows the results of Monte Carlo simulations for a homogeneous medium; 
Fig. 4 a differential backscatter signal normalized at zero absorption for several 
scattering coefficients; 

Fig. 5 shows a differential backscatter signal of a dilute suspension of 0.2 jxm 
polystyrene spheres 

30 Fig. 6 shows the total differential backscatter signal as a function of the reflection 
coefiBcient fJsO^ in the range of 10-100 mm"^; 

Fig. 7 shows the measured and calculated average path length r as a function of the 
average scattering coefficient; 



Fig. 8 is a graph of measurement results for three different absorption coefficients /4x 
showing the dc-fiber signal I, Ihe c-fiber signalJand the differential backscatter signal 
Rbs as a function of wavelength; 

Fig. 9 shows a typical* spectrum of an absorption curve A along with the specific 
S absorption coefficient of Evans Blue dye; 

Fig. 10 shows a measured -4 * as a function of the absorption coefficient /4, at 600 
nm; > 

Fig. 1 1 shows typical spectra measured in a suspension of 1 .0 |xm polystyrene spheres 
with and without Evans Blue dye; 
10 Fig. 12 graphically shows a molar extinction coefficient as a function of wavelength; 
Fig, 13 shows in vivo measurements and a fit of the differential backscatter signal Rhs 
in a human trachea reaUzed using a fiber diameter of 400 /xm, and 
Fig. 14 shows in vivo measurements and a fit of the differential backscatter signal Rbs 
in a human trachea showing very low oxygenation indicative for lung tumor. 

15 

A schematic diagram of a preferred embodiment according to the invention is 
shown in figure 1. The setup consists of a set of optical fibers for the delivery and 
collection of light to and from a sample 1 under investigation. Light firom a light soxirce 
2, for example a Tungsten Halogen lamp (Avantes HL-2000-FHSA), is led through a 

20 first arm 3 of a bifurcated optical fiber. The bifurcated optical fiber is at a distal end 4 
coupled to a first distal end of a delivery- and-coUection fiber 5 (in the following 
referred to as dc-fiber 5) which is small enough to be fit through a working channel of a 
clinical endoscope, not shown. A second distal end of the delivery- and-coUection fiber 
5 contacts the sample 1. Alongside the dc-fiber 5, a collection fiber 6 is arranged to 

25 collect light reflected by the sample 1 . The collection fiber 6 (referred to as c-fiber 6) is 
connected to a slave channel of a dual-channel spectrometer 7, for example an Avantes 
SD2000. Preferably, the dc-fiber 5 is polished at a small angle to reduce specular 
reflections... ... 

Light reflected back from the sample 1 into the c-fiber 6 is led direcfly into the 

30 slave chaimel of the dual-channel spectrometer 7. A second arm 8 of the bifiircated 
fiber is connected to a master chaimel of the dual-chaimel spectrometer 7. Light 
reflected into the dc-fiber 5 is coupled back into the bifurcated fiber, and reaches the 
dual-chaimel spectrometer 7 via the second arm 8 of the bifurcated fiber. An output of 



the spectrometer 7 is connected to an input of a processor 9 which is arranged to 
analyze signals ficom fhe spectrometer 7. 

If only fhe dc-fiber 5 would be used to deliver and collect light to and from the 
sample 1, a large fraction of collected light is due to single backscattering from small 

5 sample depths, see [1]. A single to multiple scattering ratio depends on the scattering 
coefficient and phase function of the sample 1 and on a diameter of the dc-fiber 5. The 
contribution of multiply scattered light to the signal of the dc*-fiber 5 can be 
approximately determined by combining the signal of tiie dc-fiber 5 with a signal 
coming from an additional fiber, i.e. the c-fiber 6 mentioned above. 

10 In [4] a differential backscatter signal jR^^ as a function of the wavelength X is 

determined using a formula like 



J 



(I) 
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where IQ^ is the signal from the dc-fiber S in contact with the sanople 1, In(}) is the 
signal from the dc-fiber S submersed in a fluid with an appropriate iefi:active index (for 
tissue: water would be appropriate), ly^hue 09 is the signal from the dc-fiber 5 with fhe 
probe-tip at a specific distance fix>m a diffuse reflecting reference material with a large, 

20 preferably wavelength-independent reflectance coefficient (white spectralon) and Ibiadt 
(>) is the signal &om tiie dc-fiber 5 with the probe-tip at that same specific distance 
from a diffuse reflecting reference material with a small, preferably wavelength- 
independent reflectance coefficient (black spectralon). Furthermore, J ^ is the signal 
from the c-fiber 6 in contact with fhe sample 1 and Jwhite/biaek 09 is the signal fix>m fhe c- 

25 fiber 6 witii the probe-tip at the previously mentioned specific distance from tiie 

white/black spectralon. Finally, c is a calibration constant that depends on the distance 
.between the probe-tip and the reference .material^. . , 

According to tiie invention, the processor 9 is arranged to calculate the physical 
feature using a predefined math^atical model, fhe differential backscatter signal (Rbs) 

30 and a curve fittmg mechanism. In an embodiment, the diameter of the fibers S, 6 are 
selected depending on a mean free path (raQ)) of photons sent into fhe sample 1 . It is 
noted that if the mean free path can not be estimated before selecting a fiber diameter. 



initially two arbitrary fiber diametCTS may be selected. After curve fitdng the measuring 
results using two different mathraiatical models, it will show which model applies. 

Figures 2a and 2b depict fiber tips of the dc-fiber 5 and the c-fiber 6 in the 
situation wherein the mean firee palh (mfp) of photons con:iing out of the dc-fib^ 5, is 

5 much larger than a diameter dfiher of the fibers 5, 6. In an embodiment, flie diameters of 
both fibers 5, 6 are of equal size, however it should be understood that other selections 
are possible. In figure 2a, lines 21 and 22 show an example of a path traveled by a 
detectable singly scattered photon. In figure 2b, liaes 23, 24, 25 and lines 23, 24, 26 
show two possible paths of detectable multiply scattered photons. All multiple 

10 scattering events occur at such large -distances £com the fiber tip of the fibers 5, 6, that 
the probability of detection of multiply scattered photons is roughly eqtial for the dc- 
fiber 5 and the c-fiber 6. The differentia] backscatter signal RbsW will now purely be 
determined by singly scattered photons. 

15 In an embodiment, the diameter of the fibers 5, 6 are selected so that m£p>dfiber' 

In the predefined mathematical model of this embodiment, the differential backscatter 
signal Rbs(X) is an exponential function of two times the mean free path. Below, an 
explanation for this model is given. 

In the absence of absorbers, the differential backscatter signal RbsO^ is 

20 proportional to the local, superficial scattering coejBBcient jLhO^'^QscaW p * As : 

i?^.rV=Capp~- |dn./^(A,n)-G,^(A)-p-As (2) 

where Capp is an apparatus constant that depends amongst others on the distance 
25 between the probe tip and the reference materials (black and white spectralon), 
is a function called the phase function where Q is the scattering angle, QscaO^ the 
scattering efficiency, p fee ooncelntration of substaiices present in thfe sample 1, arid As 
tiie area of a scattering particle. For example, using a fused silica fiber with numerical 
aperture NA^^ 0.22, the differential backscatter signal RbsO^ can be approximated by 



30 



RbsO^" Capp • f- Jd(p - Jd9 . sm(e). p{WO) ■ fi, (A) = C^'- im^sO^ (3) 



4^ 0 no 
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where <j} is the azimuthal angle and 9 is the polar angle. 

Figure 3 shows the differential backscatter signal ii^^^T^ of measurements (see 
5 dots) of a dilute suspension of 0.2 jiim polystyrene spheres along wifli a calculation (see 
curve 32) according to Eq. (3). In figure 3 RbsO^ is shown using arbitrary units (a.u.). 
Also a value of Qradar where Qradar 47C- i?(\ 1 SOyOscaO^, is indicated in the figure. 
Figure 3 shows excellent agreement between the measurement (i.e. the dots) and the 
calculation, which indicates that if mfp>d/iber , the single scattering is indeed the 
10 dominant contributor to the differential backscatter signal JR^sO^ as defined in Eq. (1). 

A singly scattered photon first travels firom the tip of the dc-fiber 5 to a particle, 
and then (the same distance) back from the particle to the tip of the dc-fiber 5 (or tip of 
c-fiber 6), see also figure 2a. So an average path length r()^ traveled by the measured 
single scattered photons is equal to two times the mean free path mj^(X), i.e. 

15 

^2'mfp(h) (4) 

In the presence of « absorbing species with specific absorption coefficients jUa''''^0^, 
the dififeiential backscatter signal becomes 

20 

RbsO^ = Capp' • pQ^ 180) • /4(X) • exp(- <X) • J] pi • /fc«^'(X)) 

= Capp' • J3(\ 180) • /4(X) • exp(-2 • nift)(X)) • J] pi • fkf^^'Q^) (5) 

* 

where Capp' is an apparatus constant, p(X 1 80) is the phase function, /4(X) is the 
25 scattering coefficient of the medium, X is the wavelength of the first and second 

backscattered radiation, mj^(X) is the mean firee path as a function of the wavelength, n 
is the number of substances in the sample 1, Pi is the concentration of absorber i present 
in a detection volume of the sample 1, and fia^^^'^'^Qs) is the absorption coefficient of 
absorber z as a fimction of the wavelength. 
30 It is noted that in Eq. (5) the assumption is made that absorbers are 

homogeneously distributed and do not influence each other. The Eq. (5) may be 
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corrected for non-linear phenomena such as an inhomogeneous distribution of 
absorbers, see e.g. [8]. 

According to an embodiment, the specific absorption coefficients of the 
absorbers, the wavelength dependency of tiie scattering coefficient jUg and the phase 
5 function /?, together with Eq. (5) are used in order to calculate tiie concentrations of all 
the absorbing substances present in the detection volume of the sample 1. Since the 
detection volume is typically very small in the present invention, the extracted 
. concentrations are highly spatially resolved. This is not possible with the known 
methods that are based on diffuse reflectance, and wherein the obtained concentrations 
10 are averages over large sample volumes, see e.g. [2]. 

The apparatus constant Capp* (Eq. 3) can be determined for a specific distance 
between the tip of the dc-fiber 5 and tiie reference materials (black and white 
spectralon). For a suspension of monodisperse polystyrene spheres of known size and 
concentration, the scattering coefficient fjg and the phase fiinction /?(180) can be 
15 calculated using Mie theory [4]. The apparatus constant Capp' simply follows ftom Eq. 
(3). In terms of the volume fraction / of the suspension, the radius of the spheres a and 
the radar efficiency coefficient QradarQ^^ 47t;' pQ< ,1 ZOyQscaQ^ the apparatus constant is 
determined by 



20 i?fo(X) » Capp*- J3(\180)74y(X) = Capp ♦ 0.05968 QradatQ<) (6) 

a 



According to another embodiment, the selected diameter djwer is chosen so that 
the mean fiee path is smaller than dfiter- In this embodiment, the differential backscatter 
signal Rbs is a function of the fiber diameter dfiber* This will be discussed in more detail 
25 below. 

When the mean fi:ee path of the photons is smaller than the selected fiber 
. diameter (i.e. rsxSg(X)<dfib&)y the contribution of multiply scattered light to the 
differential backscatter signal RbsW of tihie single dc-fiber 5 caimot completely be 
removed using Eq. (1). In this case, it appears that the average path length of the 
30 photons contributing to the signal Rbs(X) becomes nearly independent of the optical 
properties of the sample 1 . In this situation, multiple scattering events already occur at 
small distances from the tip of the dc-fiber 5. An analytical expression for the 
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backscatter signal Rbs(^) is not available for fhis situation and Monte Carlo simulations 
were used to model the behaviour otRbsW as a function of tibie diameter of the fibers 5, 
6 and of the optical properties of flie sample 1. Figure 4 shows the results of Monte 
Carlo simulations using the MCML-code ^onte Carlo for Multi-Layered media) of 
5 Wang et al [6,7] for a homogeneous medium with an anisotropy value g=0.9. A flat 
circular incident beam with diameter dfiter is directed onto the sample 1, and the 
differential backscatter signal Rbs is calculated by subtracting the total reflectance in the 
c-fiber 6 (with diamet^ dfiber and center located at a distance djiber fcom the center of the 
incident beam) from flie total reflectance in the dc-fiber 5 (with diameter dfiter 
10 overlapping the incident beam). Simulations were performed for sets of four different 
scattering coefficients (|J«=15, 25, 50 and SO mm'^), four different fiber diameters 
(^^6c;^200, 400, 600 and 800 [ixn) and five different absorption coefficients (^=0, 0.2, 
0.4, 0.6 and 0.8 mm'^). 

Figure 4 shows Rbs as a fiinction of absorption coefficient ^ where the open 
IS circles/dashed lines correspond to 4^^=200 pm, the filled circles/dotted lines 

correspond to dfibot^ao |xm, the open squares/solid lines correspond to djUfei^&iQ pm 
and the filled squares/dash-dotted lines correspond to ^er^SOO yim. The differential 
backscatter signal Rbs for each scattering coefficient was normalized to the 
idfiber^QO Mm, pa'^O mm"*) case. Figure 4 shows liiat in the absence of absorption, i.e. 
20 Pa = 0, the differential backscatter signal Rbs depends linearly on the scattering 
coefficient ps. Furthermore, the slope of the strai^t lines (signifying the relation 
between Rbs and pa) depends only on the fiber diameter and is independent of the 
scattering coefficient ps. The latter is more clearly demonstrated in figure 5, where the 
differential backscatter signal Rbs is normalized to unity at zero absorption for all 
25 scattering coefficients p«. The open circles correspond to djiber'=<^^0 pm and the filled 
squares correspond to dfib&^iQO pm. These Monte Carlo simulations therefore suggest 
that inthe.situation where mQ)<dfiber» the .diameter of fiie fibers 5, 6 determines the 
average path length x of the measured photons. The backscatter signal Rbs for fhis range 
of parameters can thus be written as 

30 

Rbs(^) = Cr Ps • exp(-T • Pa) = Cr Ms • exp(-C2 • djiber • Ma) (7) 
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where Ci and Cz are constants, x is the average path lengtti, is the absorption 
coefiBcien^ liia is the scattOTng coefBcient and dju,er is the fiber diameter of the fibers 5, 
6. 

5 An exact analytical expression for RtsM is not available due to the large 

contribution of multiple scattering events to the signal. Measurements were done for 
determming a total integrated backscatter signal Rtoi for a range of A between 400-900 
nm, using the formula 

900nm 

Figure 6 shows that the integrated total backscatter signal RtotCAk) is proportional to 
fisO^ in the relevant range of /4 betweea 10-100 nmi'^ . Therefore, in the absence of 

■ 

absorbers, it follows that 

15 

Rl^Os) = CappMsO^ 

which is in agreement with the Monte Cado simulations. 

20 In the presence ofn absorbing substances in a suspension, with specific 

absorption coefficients jj/^^'^W* the differential backscatter signal becomes 

* 

RisO^ = C„^p /4fV exp(-T- i: Pi • Ma'^^W) (10) 



25 



11 

Figure 7 shows tibie measured and calculated average path length r as a function 
of the average scattering coeflacient <fisO^> (with 500 nni< X< 700 nm) for dfiher = 0.4 
mm and for absorption coefficient ~ 2,0 nrai'^ at X=600 mn. In figure 7 
measurement are indicated by dots, the r — 2*mi^ curve is indicated by a line 71 and 
Monte Carlo simulations are depicted by dashed lines. Identical results were obtained 
for suspensions with an absorption coefficient ot/ja—lS> mxxL^ at 600 nm. The average 
path length r was determined using suspensions of polystyrene sph^s with different 
sizes and concentrations to vary the scattering coefficient ids(h). The anisotropy g of 
these suspensions was in the range of 0.8-0.9. Evans Blue dye was added as an 
absorber, and the average path length r was calculated j&om Eqs. (9) and (10) and 
knowledge of tihte concentrations and specific absorption coefficient of Evans Blue, as 
will be known to the skilled person. 

Looking at the measured average path lengths of figure 7, it clearly shows that for 
large scattering coefiicients (jJk - 10-100 nun ^ the range relevant for tissue) the 
average path length r is independent of the scattering coefficient to within 10% and 
approximately equal to half the fiber diameter (r « 0.24 mm while dfiber ~ 0.40 mm), in 
agreement with the Monte Carlo simulations (the dashed lines correspond to Monte 
Carlo calculations for djiber 0.2, 0.4, 0,6 and 0.8 mm). For snaall scattering coefficients 
(e.g. /(y < 5 nrni"^) the average path length r is well described by r== 2-mfp, according 
to Eq.(4), see line 71. Figure 7 also clearly demonstrates that the ttaiosition firom the 
'single scattering regime' to the 'constant path length regime* occurs for mean free 
paths of the order of the fiber diameter. It is therefore expected that single scattering 
will prevail over a larger range of scattering coefficients for fiber diameters smaller 
than 400 pm. 

In the following, the effect of absorption on the average path length t will been 
examined in more detail. Various concentrations of Evans Blue dye were added to a 
. suspension of polystyrene spheres with scattering coefficients fJs of 35 mm**^ The 
concentrations of Evans Blue (EB) dye were varied such that the absorption coefficient 
/ja at 600 nm was in the range of 0 to 2 mm"^ Typical results of the differential 
backscatter sigoal Rta for three different absorption coefGcients fJa are shown in figure 
8. Note that the signal / of the dc-fiber 5 is plotted on a different vertical scale than the 
signal J of the c-fiber 6 and the differential backscatter signal Rts^ 
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The spectra with Evans Blue present in Hie suspension were divided by the 
spectrum with no Evans Blue present and the negative natural logarithm of the tatio 
R^/R^ was determined: 

5 ^ = -ln(R^/R^)-T-p-/A,'^^'™ (11) 

where p is the concentration of Evans Blue, and ^4,^^^^'^^ is the specific absorption 
coefficient of Evans Blue. 

10 Figure 9 shows a typical spectrum of an absorption curve 92, along with the 

specific absorption coefficient fjJP^"^^ of Evans Blue dye, see curve 94. 

For all concentrations, an area -4* under the absorption curve 92 was determined 
in the wavelength range X between 500 and 650 nm. From Eq, (1 1) it follows that if the 
average path length x is independent of the absorption coefficient jj^t^^^ area A* should 
15 depend Imearly on the concentration p of the sfpecies in the suspensioiL 

Figure 10 shows a measured * as a function of the absorption coefficient fJa at 
600 nm. A curve fitting is done, resulting in a line 104 for /4(X>^35 mm"^ Figure 10 
shows that the average path length r is indeed independent of the absorption coefficient 

fia in the range 0-2 mm'^ . 

20 From the previous results of figure 4 to 10, it shows that for m^ < rf/r6cr» the 

differential backscatter signal R^s is described by Eq. (7) with Cz « 0.6- Figure 1 1 
shows typical spectra measured in a suspension of LO pm polystyrene spheres with and 
without Evans Blue dye (|ja= 2 and 0 mm~^ at 600 nm^ respectively). From Eq. (7) and 
figure 7, the relation between the differential backscatter signals with and without 

25 absorber is given by 

Rf^(^^fJ^ = Rts(J^0);exp(-0.24'/j^ (12) 

The calculated spectrum according to Eq. (12) is plotted as a dashed line 1 10 in 
30 figure 1 1 and shows excellent agreement with the measured Rbsi^J^)* see line 1 1 1 . In 
figure 11, line 112 depicts Rbsi^O). 
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In short, the average path length x of photons measured when subtracting the 
signals of the c-fiber 6 from the dc-fiber 5 using Eq. (1), is independent of the optical 
properties of the sample 1 and approximately equal to half the diameter of the fibers 5, 
6 used, as long as the fiber diameter is larger than the m^. 



In a specific embodiment^ the device according to the invention is arranged to 
determine concentrations of oxygenated blood in tissue. Since the scattering coefficient 
of tissue ^^^^"^ is in the range of 10-100 mm"^ the fiber diameter should be smaller 
than a certain maximum diameter dmax where dmax is between 10 and 100 (mi, so for 
example smaller than SO nm, in order to measure predominantly single scattering in 
tissue. In this case, Eq. (5) holds. For fibers 5, 6 with much larger diameters (e.g. 200 
or 400 pm), the differential backscatter signal RbsQ^) is described by Eq. (10) with r » 

It is presently known that the wavelength dependence of the scattering coefficient 
in tissue /LiJ^^^ can be adequately described by an empirical power-law function, see 
also [3], [4], [5]. 

/4''^^\X)= ei'X^ (13) 

with a and b constants that depend on the size, concentration and relative refiiactive 
index of the scatterers (i.6. substances) preset in the detection volume. 
The donoinant absorbers in tissue in the visible wavelength range are oxygenated and 
deoxygenated blood. Thus in tissue Eq. (10) becomes 

RM = Capp • aX^ • exp(- 0.6 dju,er • Pbiood • (S02 - /A,^"^^ + (I-S02) • )) 

= C'app • V' • exp(-. 0.6 dju^ • pblood • (S02 • M/^'^ + (I-S02) • fia'^*'^ )) (14) 

where pbiood is the concentration of blood, S02 is the blood oxygenation (percentage 
oxygen saturation) in a certain detection volume, Capp is a constant that depends on the 
caUbration constant c, Capp is Capp-a, X is the wavelength, b is the slope of the 
scattering coefficient defined in Eq. 13, fjj^^^^ is the specific absorption coefficients of 
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folly oxygenated blood, fj^P^'^"^"" is iJie specific absorption coefficients of folly 

deoxygenated blood. 

Non-linear phenomena such as an inhomogeneous distribution of absorbers are 
not incorporated in Eq. (14), but can be added by the skilled person, see e.g. [8], 

Since the specific absorption coefficients of folly oxygenated (j^a^"^ and folly 
deoxygenated Qj/^'^''^^ blood are weU known, see figure 12, Eq. (14) can be fitted to 
the measured data to yield the slope b of the scattoing coefficient /Js^"^\ the ^ 
concentration pwoodand the oxygen saturation S02 of the blood present m the detection 
volume. When a correction is made for the inhomogOTeous distribution of blood in the 
vessels, a vessel diameter D may be determined as well. Since the average detection 
depth is small (e.g. 0.1 mm), the blood present in Ihe detection volume when measuring 
non-invasively is located in capillaries. 

In figure 13, in vivo measurements of backscattering in a human trachea together 
with a fit using Bq. (14) are shown. The measurements are realized using a fiber 
diameter of 400 jum. The dots depict Ihe measurements and curve 130 is a fittmg curve. 
In figure 13, b - -0.94 and the oxygenation S02- 95%. 

The present invention can be used for tumor detection. Tumor growth may, due 
to its excessive oxygen consumption, be accompanied by a low capillary oxygen 
saturation, which can only be assessed using a very localized measurement. Since (pre- 
)cancerous tissue is generally more heterogeneous than normal tissue, the standard 
deviation of multiple measurements is likely to be larger for (pre-)cancerous tissue than 
for normal tissue. Standard deviations in the measurements can be calculated for the 
oxygen saturation, the blood concentration, the blood vessel diameter, and the slope 6 
of the scattering coefficient /j/^^. It is noted that the invention is by no means 
restricted to determine a concentration of a substance as the physical feature. All 
features, mentioned in the previous phrase can be regarded as physical features. 

An exaicnple of a measurement of a Ixmg tumor is shown in figure 14. The- shape 
of the dip in the wavelengfli range of 500-600 nm in this figure demonstrates the 
depletion of oxygen from the capillaries of this tumor due to its excessive oxygen 
consumption. 



15 

When a needle-probe is used^ the local oxygenation and scattering coefficient 
Z^"**"^ can be measured invasively. This could be helpful in determining tumor-margins 
intra-operatively in real-time, for instance during resection of a breast-tumor. 
According to an embodiment, the device comprises multiple probes and a 
5 multichaxmel spectrometer for multiple simultaneous measurements on different 
locations of the sample 1. Using this device, multiple measurements can be made 
simultaneously on different locations of for example a suspicious lesion. 

In yet another embodiment, the device comprises at least two pairs of fibers, 
having different fiber diameters. For example, when a pair of fibers with 100 ftm, a pair 
10 of fibers with a diameter of 200 jum and a pair of fibers with a diameter of 400 /xm axe 
used, information firom different depths in the sample 1 can be obtained as the average 
pathlength increases with increasing fiber diameter. 

The method and apparatus according to the invention can also be used to analyze 
local drug concentrations. From Eq. (10) it follows that if the specific absorption 
IS coefficient of a certain dmg is known, the local concentration p of that substance can 
be detemiined using the invention. 

Another possibility of the present invention is to monitor glucose concentrations. 
The scattering coefficient pj'^^^^ depends among others on the relative refiractive index 
of the scatterers with respect to the surroimding medium (in tissue: cytoplasm). The 
20 reftactive index of the surroimding cytoplasni is likely to depend on the concentration 
of glucose. A change in the glucose concentration will therefore likely affect the slope 
. b of the scattering coefficient fij^^^ , see Eq. (13). 

Whilst specific embodiments of the invention have been described above, it will 
be appreciated that the invention may be practiced otherwise than as described. For 
25 example, a concentration of a substance in pointed water may be calculated. The 
description is not intended to limit the scope of the invention. 



30 
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Claims 

1 . A device for the detennination of a physical feature of a medium^ comprising: 

- a light source (2) for producing radiation; 

5 - a probe with at least a first and a second optical fiber (5, 6), said first optical fitber (5) 
having a first diameter and being ananged to deliver said radiation on a sample (1) of 
said medium and to collect first backscattered radiation fi:om said sample (1), said 
second optical fiber (6) having a second diameter and being airanged to collect second 
backscattered radiation, wherein said second optical fiber (6) is positioned alongside 

10 said first optical fiber (5); 

- a spectrometer (7) for producing a first signal (I) based on said first backscattered 
radiation, and for producing a second signal (J) based on said second backscattered 
radiation; 

- a processor (9) arranged to determiae a differential backscatter signal (Rbs) using said 
1 5 first and second signals (I, J), 

characterized in that said processor is arranged to calculate said physical feature using a 
predefined mathematical model, said differential backscatter signal (Rbs) and a curve 
fitting mechanism. 

20 2. Device according to claim 1, wherein said first and second diameter are selected 
depending on a mean firee path of photons sent into said sample (1). 

3. Device according to claim 1 or 2, wherein said first and second diameter are 

« 

essentially equal to a selected fiber diameter. 

25 

4. Device according to any of the preceding claims, wherein said physical feature is 
a concentration of at least one substance in said medium. 

• • 

5. Device according to claim 3 or 4, wherein in said predefined model said 

30 differential backscatter signal (Rbs) is an exponential function of two times said mean 
firee pafii, and wherein said selected diameter is selected so that. 

where dfiber is said selected diameter and mfp is said mean free path. 
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6. Device according to claim 5, wherein said differential backscatter signal (Rta) is 
defined by: 

-Rte^V = Capp' pCMSO) ■ fJ,(K) ■ exp(-2 • ni^(X)) • ^ Pi " /«.^'0s)) 

5 

where Capp' is an apparatus constant, p(\ 1 80) is a phase function, /^(X) is a scattering 
coefficient of said medium, Xis the wavelength of said first and second backscattered 
radiation, mfp(X) is said mean free path as a function of the wavelength, n is the 
niunber of substances in said medium, pi is the concentration of absorber i present in a 
10 detection volume of said sample (1), and jjJ^^^'^Q^ is an absorption coefficient of 
substance i as a function of the wavelength. 

7. Device according to claim 3 or 4, wherein in said predejQned model said 
differential backscatter signal Rbs is a function of said selected fiber diameter, and 

1 5 wherein said selected diameter is selected so that: 

dfiber>raSp 

where dja„r is said selected diameter and mJ^ is said mean free path. 

8. Device according to claim 7, wherein said differential backscatter signal Rbs is 
20 defined by: 

Rbs = Cr He ' e3qp(-C2 • • Ha) 
where Ci and C2 are constants, |ia is an absorption coefficient of a medium, |is is said 
scattering coefficient and dfiter is said selected fiber diameter. 

25 9. Device according to claim 8, wherein C2 is approximately 0.6. 

10. DeAdce according to any of the preceding claims, wherein said medium is living 
tissue. 
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1 1 . Device according to any of the preceding claims, wherein said device comprises 
multiple probes and a multichatmel spectrometer for multiple simultaneous 
measurements on different locations of said sample (1). 



S 12. Device according to any of claims 1 »1 0, wherein said probe comprises at least 
two pairs of fibers, having different fiber diameters. 

13. A method for the determination of a physical feature of a medium, comprising: 

- producing radiation with a light sowce (2); 

10 - placing a probe with at least a first and a second optical fiber (5, 6), said first optical 
fiber (S) having a first diameter, on a sample (1) of said medium; 

- placing a second optical fiber (6) having a second diameter, alongside said first 
optical fiber (5); 

- sending light coming firom said light source, througjh said first optical fiber; 

IS " collecting first backscattered radiation through said first optical fiber and second 
backscattered radiation through said second optical fiber; 

* producing a first signal Q) based on said first backscattered radiation, and a second 
signal (J) based on said second backscattered radiatioi^ 

- determine a differential backscatter signal (Rbs) using said first and second signals (I, 
20 J), 

characterized bv 

- calculating said physical feature using a predefined model, said differential 
backscatter signal (Rbs) and a curve fitting mechanism. 



25 14. Method according to claim 13, wherein said first and second diameter are 
selected depending on a mean free path of photons sent into said sample (1). 



IS. Method according to claim 13-14, wherein said method comprises: 

- simultaneously measuring backscatter radiation on different locations of said sample 
30 (1); 

- determining a physical feature for said different locations; 

- calculating a standard deviation of said physical feature. 
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1 6. Computer program product to be loaded by a computer, said computer program 
product, after being loaded, providing said computer with the capacity to: 

- receive a first signal (I) indicative for a collected radiation received from a first fiber 
(5) and a second signal (J) indicative for a collected radiation received from a second 

5 fiber (6); 

- receive a selected fiber diameter parameter indicative for a diameter of said first and 

second fiber (S, 6); ^ 

- deteraiine a differential backscatter signal (Rbs) using said first and second signals Qi, 

J); 

10 - calculate a physical feature using a predefined model, said differential backscatter 
signal OElte), a curve fittmg mechanism and said selected fiber diameter parameter, 

1 7. Data carrier provided with a computer program product according to claim 1 6. 

15 



1 
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Abstract 

The invention relates to a device for the determination of a physical feature of a 
medium^ such as a concentration of a substance in the medium. The device comprises a 

5 a probe with at least a first and a second optical fiber. The fixst fiber is used to deliver 
and collect radiation to and firom the sample. The second fiber is used to only collect 
radiation. Radiation collected by the two fibers is transferred by a spectrometer that is 
connected to a procesor. Depending on the ratio between the diameter of the fibers and 
a mean firee path of photons sent into the sample^ two different mathematical models 

10 are used by a program to be loaded by the processor. The models show very good 

agreement with in vivo measurements. With the device, very accurate measurements of 
the concentration of substances in the sample can be done. In addition, these 
measurement can be done very locally because the two fibers are positioned alongside. 

15 [Figure 1] 
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